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1. Introduction 

Understanding the behavior of large pools of credit assets is currently a problem of central importance. 
Banks often hold such large pools and their risk-reward characteristics need to be carefully managed. In 
many cases, the losses in the pool are (hopefully) rare as a consequence of diversification. In |Sowa| . we have 
used the theory of large deviations to gain some insight into several aspects of rare losses in pools of credit 
assets. Our interest here is the effect of recovery. While a creditor either defaults or doesn't (a Bernoulli 
random variable), the amount recovered may in fact take a continuum of values. Although many models 
assume that recovery rate is constant — i.e., a fixed deterministic percentage of the par value, in reality the 
statistics of the amount recovered should be a bit more complicated. The statistics of the recovered amount 
should depend on the number of defaults; a large number of defaults corresponds to a bear market, in which 
case it is more difiicult to liquidate the assets of the creditors. Our goal is to understand how to include this 
effect in the study of rare events in large pools. We would like to look at these rare events via some ideas 
from statistical mechanics, or more accurately the theory of large deviations. Large deviations formalizes 
the idea that nature prefers "minimum energy" configurations when rare events occur. We would like to see 
how these ideas can be used in studying the interplay between default rate and recovery rate. 

Our work is motivated by the general challenge of understanding the effects of nonlinear interactions 
between various parts of complicated financial systems. One of the strengths of the theory of large deviations 
is exactly that it allows one to focus on propagation of rare events in networks. Our interest here is to see 
how this can be implemented in a model for recovery rates which depend on the default rate. 

This work is part of a growing body of literature which applies the theory of large deviations to problems of 
rare losses in credit assets and tries to approximate the tail distribution of total losses in a portfolio consisting 
of many positions. For a survey on some existing large deviations methods applied to finance and credit 
risk see |Pha07| . In |GKS071 IGL05| rare event asymptotics for the loss distribution in the Gaussian copula 
model of portfolio credit risk and related importance sampling questions are studied. In |G02j the author 
considers saddle point approximations to the tails of the loss distribution. Measures such as Value-at-Risk 
(VaR) and expected shortfall have been developed in order to characterize the risk of a portfolio as a whole, 
see for example [DS03| [G03l IMFE051 IFM02[ IFM03| . Issues of recovery have also been considered in the 
works [XBRS05 , SH09j and in the references therein. Our work is most closely related to |DDD041lLMS09j . 
where the dynamics of a configuration of defaults was studied. In |DDD04j it is assumed that a "macro- 
environmental" variable Y, to which all the obligors react, can be chosen so that conditional on Y, the 
recovery and default rates in a pool of finite number of "types" are independent. In jLMSOQj the recovery 
rate is assumed to be independent of the defaults in a pool of one type with a brief treatment of a pool with 
finite types. Our work here is explicitly interested in the dependence of the recovery rate on the fraction of 
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the number of defaults and the framework of our efforts is a continuum of types. The case of a continuum 
of types requires shghtly more topological sophistication. The dependence assumption and the continuum 
of types unavoidably complicate the proof since several auxiliary technical results need to be proven. 

The novelty of our result stems from the following two things. The first novelty is that the distribution 
of the recoveries for the defaulted assets depends upon the number of defaults in a fairly general way. This 
allows for consideration of the case when recovery rates arc affected by the number of defaults. Note however 
that the individual defaults are assumed independent. The second novelty is more technical in nature. In 
particular, we prove the large deviations principle for the joint family of random variables (Ljv , vn ) where L n 
is the average loss and i^n is the empirical measure on type-space determined by the names which default. In 
addition to the large deviations principle, we derive various equivalent representations for the rate function 
which give some more insight to the favored way to rearrange recoveries and losses among the different types 
and also ease its numerical computation. 

The paper is organized as follows. In Section [2] we introduce our model and establish some notation. 
In Section |3] we study the "typical" behavior of the loss of our pool; we need to understand this before 
we can identify what behavior is "atypical" . In Section |31 wc present some formal sample calculations and 
our main results. Theorems 14.61 and 14.81 These calculations are indicative of the range of possibilities and 
they illustrate the main results. In Section [S] we present some numerical examples which illustrate some of 
the main aspects of our analysis and conclude with a discussion and remarks on future work. The proof 
of Theorem 14.61 is in Sections [51 [7] and |S1 Section [S] contains an alternative expression of the rate function, 
which is a variational formula which optimizes over all possible configuration of recoveries and defaults, and 
which leads to a Lagrange multiplier approach which can be numerically implemented. Lastly, Section [lU] 
contains the proofs of several necessary technical results. 

The model at the heart of our analysis is in fact very stylized. Since our primary interest is the interaction 
between default rates and recovery rates, our model focuses on this effect, but simplifies a number of other 
effects. In particular, we assume that the defaults themselves are independent. Our work complements the 
existing literature and hopefully contributes to our understanding of the interplay and interaction of recovery 
and default rates. It seems plausible that more realistic models (e.g., which include a systemic source of 
risk) can be analyzed by techniques which are extensions of those developed here. 

2. The Model 

In this section we introduce our model and explain our goals. Let's start by considering a single bond (or 
"name"). For reference, let's assume that all bonds have par value of $1. If the bond defaults, the assets 
underlying the bond are auctioned off and the bondholder recovers r dollars, where r G [0, 1]. We will record 
the default/survival coordinate as an element of {0, 1}, where 1 corresponds to a default and to survival. 

The pool suffers a loss when a bond defaults, and the amount of the loss is $1 — r, where r £ [0, 1] is the 
recovery amount (in dollars). Define 

^ dcf ( 1 if the ri*'' name has defaulted 
" 1 otherwise 

and let f„ be the non-zero loss that occurs in the event of default of the n**^ name, i.e., when A„ = 1. 
Equivalently, £n is the non-zero exposure of the n'^ name that occurs when A„ = 1. Clearly if r„ is the 
recovery of the n*^ name then = 1 — r„. 

The default and loss rates in the pool are then 

N N 

(2.1) ^^'=^'^E^" L^<i=^fi-^£„A„. 

n— 1 n— 1 

We denote by p^'" the risk-neutral probability of default of the n"^ name in a pool of N names. Moreover, 
let gi^^"(£'Af, dr) be the distribution of the recovery of the n"^ name in a pool of names. Notice that we 
allow p^'" (ZJjv , dr) to depend on the default rate D^. 

Let us make now the aforementioned discussion rigorous by introducing the probability space that we 

consider. The minimal state space for a single bond is the set Eo {0} U ({1} x [0, 1]). Since we want to 
consider a pool of bonds, the state space in our model will be E i?^, where N =^ {1, 2 . . . }. E'o, as a closed 
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subset of {0, 1} X [0, 1] is also Polish, and thus E is also Polish. We endow E with the natural cr-algebra 
J^=^MiE). 

The only remaining thing to specify is a probability measure on (E,^). For each TV e N, fix {p^'" : 
n E {1,2...A^}} C [0,1]. These are the risk-neutral default probabilities of the names when the pool 
has N names. We next fix {g^'" : n e {1,2... TV}} C C([0, 1]; ^[0, 1] j]; i.e., a collection of probability 
measures on [0,1], the range of the recovery, indexed by the default rate Djy. Namely, for each N E N 
and n G {1, 2, • • • , N}, is a probability measure that depends on the default rate D^- Some concrete 
examples are in Subsection 15. II 

For each n G N and u! ~ {uji,lj2 ■ ■ ■), define the coordinate random variable ~ aj„. For each N G N, 

we then fix our risk- neutral probability measure P^v G 3^{E) (with associated expectation operator Ejv) by 
requiring that 



N 



n UXn) 



= E 



N 



n I (1 - P^'nUQ) + P^^" / /„(1, r)e^^"p^, dr) 
1 I Jre[04] 



for all {/n}„=i C B{Eo). In particular, the aforementioned construction implies that for each n G N, 
{A„}^^]^ is an independent collection of random variables with W'm{^n = 1} = p^'" under Pjv. 

With this probabilistic structure in place, we will clearly want to be able to condition on the default rate 
so that we can then focus on the recovery rates. To this end we define the a-algebra 

^ =V{A„ : n e N}. 

We are interested in the typical behavior and rare events in this system. We seek to understand the 
structure of these rare events in our model. We are interested in the following questions: 

- dcf 

• What is the typical behavior of the system? In other words we want to characterize L = Yurn^oo Lm ■ 

• What can we say about rare events in this general setting? In particular we compute the asymptotics 
of PatILat > /} as oo, particularly for I > L. Then {i„ > 1} is an "atypical" event. 

Understanding the structure of rare events may give guidance in how to optimally control how rare events 
propagate in large financial networks. 

Remark 2.1. Our formulation allows for a fairly general dependence of the recovery distribution on the 
default rate. This is also partially exploited in Subsection I5.il where we study several concrete examples. 
On the other hand, the general treatment of this paper unavoidably complicates the proof of the law of large 
numbers and especially of the large deviations principle. 

3. Typical Events: A Law of Large Numbers 

Let's start our analysis by identifying the "typical" behavior of Lat as — > oo defined by equation 
(|2.ip . In order to motivate the calculations, we first investigate how and behave as A'^ — >■ (X) in a 
homogeneous pool, Example 13.11 and in a heterogeneous pool of two types. Example 13.21 

Example 3.1 (Homogeneous Example). Fix p E [0,1] and g E C([0, 1]; ^[0, 1]) and let p^'^ = p and 
g^'" = g for all N eN and n e {1,2... TV}. 

Let's see what Dm and Ln look like in this setting. By the law of large numbers limjv^oo Dn = p. Thus 
in distribution 

l<n<pN 

where the ^„ 's are i.i.d. with distribution g(p, •). We should consequently have that 

lim Ln — p / {1 ~ r)g{p,dr). 



L 



N 



rG[0,l] 

□ 



^We shall write g e C{[0, 1]; .iP[0, 1]) as a map from [0, 1] X ,«[0, 1] into [0, 1] such that for each D G [0, 1], A q(D, A) is 
an element of £*^[0, 1] and such that D c>{D, ■) is weakly continuous. 
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Example 3.2 (Heterogeneous Example). Fix pA and ps in [0,1] and fix qa and qb in C([0, 1]; ,^[0, 1]). 
Every third bond will he of type A and have default probability pA and recovery distribution governed by qa, 
and the remaining bonds will have default probability ps and recovery distribution governed by qb- In other 
words, the rate of default among the A bonds is pA and the rate of default among the B bonds is pB- 
For future reference, let's separate the defaults and into the two types. If we define 

l<n<N L ' ^ 1<71<N 

neSN n^3N 

we have 

Thus liniAr^oo Dn = Pa and linijv-i.oo Dn =^ Pb, so 

defPA , 2pB 



lim Dn ^ D 

N^oo 3 3 



Thus we should roughly have 



[l<ri<p^Ar/3 l<n<2pBN/3 

where the 's have law qa{D, ■), the 's have distribution qb{D, ■) and the 's and 's are all indepen- 
dent. Consequently, it seems natural that 

\imLN = ^[ {l-r)QA{D,dr) + ^ [ {1 ~ r)gB{D,dr), 

N^co S JrelQ,l] Jre[0,l] 

the first term being the limit of the losses from the type A names, and the second term being the limit of the 
losses from the type B names. 

□ 

In view of our examples, it seems reasonable that we should be able to describe the average loss in the 
pool in terms of a frequency count of 

We note that 7^^" takes values in the set X [0,1] x C([0, 1]; ^[0, 1]). Since ^[0,1] is Polish, so is 
C([0, 1]; ^[0, 1]), and thus X is Polish. We will henceforth refer to elements of X as types. For each iV G N, 
we now define Uat G 0^{X) as 

1 ^ 

Assumption 13.31 below is our main standing assumption. It is a natural assumption since it makes the 
family {UArjArgN relatively compact. It will be assumed throughout the paper, even though this may not be 
always stated explicitly. 

def 

Assumption 3.3. We assume that U ~ limjv^oo Uat exists {in ^(X)). 

□ 

Remark 3.4. In the case of Example \3.1\ we have that U = while in the case of Example \ 3.2\ we 

have that U = |^(p^,e^) + jS^pB^es)- 

We can now identify the limiting behavior of L^. Define 

(3.1) D=^f pU{dj) and L=^f \pf {1 ~ r)g{D, dr)\ \J{dj). 



h={p,Q)ex J'(=(p,e)ex [ Jre[o,i] 

To see that these quantities are well-defined, note that 

(3.2) {p, g) ^ p and ly ^ (1 - r)v{dr) 
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are continuous mappings from, respectively, X and ^[0, 1], to [0, 1] C M. The continuity of the first map of 
(|3.2[) implies that D is well-defined. Combining the continuity of both maps of (|3.2p . we get that the map 



{p,q) (1 -r)g{D,dr) 

Jre[0.1] 

is a continuous map from X to [0, 1] C M; thus L is also well-defined. 

Lemma 3.5. Let Assumvtion \3.3\ hold and consider L given by Then, for each e > 0, we have that 



lim Vn\\Ln-L\ >e| =0. 



Proof. For g € C([0, 1]; ^[0, 1]) and D G [0, 1], let's first define 

dof 



□ 



TiD,g) 



re[0,l] 



il-r)g{D,dr). 



Again we can use (|3.2[) and show, by the same arguments used to show that D and L of p.ip are well-defined, 
that r is well-defined, and furthermore that it is continuous on [0, 1] x C([0, 1]; ^[0, 1]). For each N £ N, 
define 



N N 



dcf 1 

N 



n=l "'t 



7=(p,e)ex 
pTiD,g)UNid-f); 



7=(p,e)GX 



note that ijv is a random variable but and L^r are deterministic. Note also that by weak convergence, 

lim Ln = L and lim Dn = D. 

N~¥oo N^oc 



The claim will follow if we can prove that 



(3.3) 



lim 

N-^oc 



= and lim Ejv 

N^oo 



\Ln — Ln\ 



= 0. 



To see the first part of p.3p . we calculate that 



Conditioning on ^ = (t{A„ : n e N}, we have that 



1 ^ 



Ln — Ln 



re[0,l] 



((1 - r) - T{Dn, g'^ny g^'^'^iDN, dr) Al < 



To see the second part of p.3p , we write that Ln ~ Ln = E]^ + where 



-y 

n=l 

This implies the first part of p.3p . 

ri=l 
1 ^ 



TV 



We first calculate that 



n=l 



N 



E 



l/^^ I 



Thus liniAr_,ooIEAr = 0. Similarly, 



N 



E 



N 



Therefore 



\Dn-Dn\ 



lim 



iV2 



< 



n=l 



AN' 



Dn-D\ 



= 0. 



To bound fix > 0. Due to Assumption 13. 3[ Prohorov's theorem implies that {(Jn', N £ N} is tight, 
there is a Kr, CC C([0, 1]; ^[0, f]) such that 



Defining 



supUjv([0,l] X K^) < r]. 



Wr,{5)= sup \V{D^,q)^V{D2,q)\ 

Q&K„ 
Di.Dasio,!] 
|I5i-r>2|<(5 



for all 5 > 0, compactness of K^i and [0, 1] and continuity of F imply that lim^x^o '-^rii^) = 0. Thus 
\Sn\ < Eat 



1 

TiDN,g)-r{DN,g)\UN{d-/) 



n=l 



= E 



■'7=(p,e)6X 

<2rj + 2uj^{6) + Fn {\Dn -D\>5}. 

Take N /• oo, then let 5 \0. Finally let 77 \ 0. We conclude that indeed limTv-i-oo EAr[|£j^|] = 0. Hence, 
the second part of (|3.3p also holds which concludes the proof of the lemma. □ 

4. Problem Formulation and Main Results 

Let's now set up our framework for considering atypical behavior of the Lat's; i.e., large deviations. We 
motivate the main result, Theorem 14.61 through a formal discussion related to the setting of Examples 13.11 
and 13.21 The proof of Theorem l4.6l is in Sections [H [7] and HI In Section |6] we prove compactness of the level 
sets of the rate function. In Sections [7] and [8] we prove the large deviations lower and upper limit respectively. 
An equivalent insightful representation of the rate function is given in Theorem l4.8l and its proof is in Section 
[HI This representation is a variational formula which optimizes over all possible configuration of recoveries 
and defaults, and which leads to a Lagrange multiplier approach which can be numerically implemented. 

For the convenience of the reader we here recall the concept of the large deviations principle and the 
associated rate function. 

Definition 4.1. If X is a Polish space and P is a probability measure on (X, ^^{X)), we say that a collection 
($n)neN of X -valued random variables has a large deviations principle with rate function / : AT — > [0,oo] if 

def 

(1) For each s > 0, the set ^(s) — {x £ X : I{x) < s} is a compact subset of X . 

(2) For every open G C X , 



(3) For every closed F C X , 



lim — In^ 



lim — In ] 

n/'oo n 



'{in e G} > - inf I(x) 



'UneF}<- inf I{x). 



One thing which is clear from Example 13.21 is that we need to keep track of the type associated with each 
default (but not the types associated with names which do not default). To organize this, let ^i(X) be the 
collection of measures v on (X, ^(X)) such that i^(X) < 1 (i.e., the collection of subprobability measures). 
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As it is discussed in Sectional ^i(X) is a Polish space. For each N gN, define a random element i/^ of 
^i(X) as 

1 ^ 

(4.1) l^N ^ 5] ^"'^7"- 

n=l 

Thus, is the empirical measure on type-space determined by the names which default. 
Next, define a sequence (ZAr)jvGN of [0, 1] x ^i(X)-valued random variables as 

Since [0, 1] and ^i(X) are both Polish, [0, 1] x^i(X) is also Polish. We seek a large deviations principle for 
the Zn's. Since projection maps are continuous, the contraction principle will then imply a large deviations 
principle for the Lat's. Note furthermore that 

Dn = t^NW; 

the map v i— > i^(X) is continuous in the topology on ^i(X), so the recovery statistics depend continuously 
on vn. 

Before proceeding with the analysis, we need to define some more quantities. 



Definition 4.2. For p e [0,1], defi 



ne 



, , . def 
np{x) = 



x In ^ + (1 — x) In for x and p in (0, 1) 

Ini /or X = 1, (0,1] 

In^ for x = 0,p e [0,1) 

oo else. 



For V G .^i(X), define 

H{.) /A=(P..)^X^P (iu(7)) U(d7) - « U 
loo. else 

□ 

To understand the behavior of L^r, we need to construct its moment generating function. For this purpose, 
we have the following definition. 

Definition 4.3. For v G .^i(X), define 

K{e) =^ / J In / e'^^^-'-^g{v{X),dr) \ iy{d-f) deR 



7=(p,e)ex [ Jre[Q,i 

A:(£) =^sup{0^-A,(0)} ^e[0,l] 

sgk 

A*(£) is the Legendre-Fenchel transform of A,^{9). 

□ 

Let's now see if we can formally identify a large deviations principle for the Z]\j^s in the setting of Examples 
I3.1l and [3.2l In both examples we want to find a map / : [0, 1] x .^i(X) — > [0, oo] such that, at least formally. 



N 



{Zat e dz*} X exp[-7V/(z*)] A^^oo 



for each fixed z* = (£*, i^*) G [0, 1] x ^i(X). Assume that z^* ^ U. This is done without loss of generality, 
since, as we shall also see in the sequel, the rate function I{z*) = oo if i^* <4; U. 
Conditioning on the event that ~ v* , we have 

(4.2) Pwi^TV e dz*]=f{LN e dt\vN ^v*]fN{vN e dv*}. 

We want to derive asymptotic expressions for the first and second term in the right hand side of (|4.2p . 



Example 4.4. [Homogeneous Examvle \3.1\ continued]. First, we understand the second term in (|4.2p . i.e., 
Pjv {vn G dv*} . Recall that 

l<ri<Af 

This is now essentially the focus of Sanov's theorem-i.i.d. Bernoulli random variables. Namely, recalling 
the definition of hpix) and H{v) by Deflnition \4.2\ we have 

(4.3) VnIi^n e diy*} x exp [-Nhp (i^*(X))] = cxp [-NH (i^*(X))] 

Second, we understand the first term in (^T^, i.e. P {Lat G di*^!^^ ~ i'*} . As in Examvle \3.2l we should 
have in law 



(4.4) 



1 

N 



l<n<Ariy*(X) 



where the ^„ 's are i.i.d. with common law g{i'{X),-). In the homogeneous pool of this example, the log 
moment generating function defined in Definition \4.S\ takes the form 



A,.(0) = z/*(X)ln 



9(l~r) 



'•£[0,1] 



g{v{X),dr) 



for all 9 G M.. The logarithmic moment generating function of Ln of (j4.4p is approximately NK^jtiO/N). 
Thus we should have that 

Pat {Ln e dt\vN w V*) ^ cxp [-7VA*.(^)] 
We should then get the large deviations principle for Zj\j by combining this, (|4.3p . and (j4.2p . i.e., 

Pn {Zn e dz*} X exp [-N[Al, {£) + H{i'*)]] . 

□ 

Example 4.5. [Heterogeneous Examvle \ 3.2\ continued]. We proceed similarly to Example \4-4\ First, we 
understand the second term in (j4.2p . i.e., ¥^{1^^ e di/*}. Observe, that we explicitly have 



Similarly 



(7Aj - 



dU^'^^ 2/3 



U{7a} iV/3 



l<n<Af 



(7s) 



U{7b} 2^/3 ^^^^^ 

n^3N 



Thus 



PAr{z^Ar £ dz^*} = Pat {vn{ia\ ~ z^*{7a}, !^a'{7s} « '^*{ib}} 

dvjq 



N 



dv* dvN , s dv* 

du du^^^)'^^^^)" du^^^) 



Again, this is essentially the focus of Sanov's theorem-i.i.d. Bernoulli random variables. Namely, recalling 
the definition of hp{x) and H{v) by Definition \4.^ we have 



^nIj^n € diy*} X exp 



(4.5) 



N , 



exp 



3 



dv* 



dU 



7=(p,e)ex 



2N , 



T^A -777 (7a) - -^hs -777(7- 



dv* 



dU 



dv* 
diJ 



(7) U(d7) 



exp[-7Vi7(i/*)] , 



Second, we understand the first term in ^JTS^, i.e. ¥{Ln G di*^!^^ « i^*} ■ If ~ v* , then there are 
about Nv*{'^a\ defaults of type A, and Nv*{'^b} defaults of type B. Thus, as in Example \ 3.2l we should 
have in law 



(4.6) ^^^^1 E E 



l<n<Ni^*{-fA} l<n<Niy*{-fB} 

where the 's are i.i.d. with common law QAii^O^), •) and the 's are i.i.d. with common law £»b(:^(X), •). 
In the heterogeneous pool of this example, the log moment generating function defined in Definition \4.3\ takes 
the form 

A..(0) = i^*{7A}ln / e«(i-'')£>AKX),dr)+z.*{7s}ln / e«(l-'')^,B(l.(X), dr) 

Jr£[0,l] Jre[0,l] 

for all G M. The logarithmic moment generating function of of (j4.6p is approximately NA^t{6/N). 
Thus we should have that 

Fn {Ln G dt\iyN ~ 1^*} exp [-iVA*.(£)] 
We should then get the large deviations principle for Zn by combining this, (j4.5p . and (|4.2p . i.e., 

Pat {Zn G dz*} X cxp [-iV[A:. {£) + H{iy*)]] . 

□ 

Examples 14.41 and 14.51 motivate our main result. 

Theorem 4.6 (Main). Let Assumption hold and recall Definitions \4.<i\ and \4.3\ Then, we have that 
(Zjv)ArgN has a large deviations principle with rate function 

I{z)=Kl{(:)+H{v) 

for all z = {ItV) G [0,1] x ./#i(X). Secondly (Ljv)weN ho.s a large deviations principle with rate function 

(4.7) /'(£)- ?f /(^''^^ 

i^e^/i(x) 

for all i G [0, 1] . 

□ 

Proof. Combine together Propositions 16.11 17. 3[ and 18.11 This gives the large deviations principle for {Zn). 
The large deviations principle for (L^)jvgN follows from the contraction principle and the continuity of the 
map i> I— > J^(X). □ 



Next we mention a useful representation formula for H{h') that is defined in Definition 14.21 Its proof will 
be given in Section [TUl Define 

X,{0)''^'\n{pe<' + l-p) 
for a\\ 6 gM. and p G [0, 1]. Then Ap and hp are convex duals; i.e., 

hp{x) = sup {Ox - Xp{9)} x G R 

^'^■^^ Xp{e) = sup {Ox - hp{x)} OeR 

and wc have: 

Lemma 4.7. Let Assumvtion \3.3\ hold. Then, we have that 



H{v) - sup J / 

06C(X) 



(4.9) H{v) - sup <; / <i^{i)v{dp) - / Ap(0(7))U(rf7) 

for all V G Ji\ (X) . 

□ 



One way to interpret Theorem 14.61 is that the rate functions I and /' give the correct way to find the 
"minimum-energy" configurations for atypicahy large losses to occur. In general, variational problems in- 
volving measures can be computationally difficult, so Section IH] addresses some computational issues. In 
particular, we find an alternate expression which takes advantage of the specific structure of our problem. 
Define 

M,{e,D) =^ In f e'"-^~'-^g{D,dr) 6* G M 

(4.10) ^"^^["'^l 

Ig{x, D) =^ sup {Ox - Mg(e, D)] X G M 

for all D £ [0, 1]. Observe that we can write 

(4.11) AM = / M,{0,,y{X))iy{d^) 

for aU V e .^i(X) and 6* e M. Define next B B(X; [0, 1]). 

Theorem 4.8. Let AssumvtionW^hold. For £ e [0, 1] and U e ^{X), set 

J'ii) \ [ ($(7)4 f^i>ij), I $(7)U(d7)) +ftp($(7))l U(d7) : 



hHp,e)&^ ^ \ "'Tex 

$(7)vI'(7)U(d7) = 



(4.12) 

/7ex 

We have that I'{i) ~ J'i^) for all £ G [0, 1]. An alternate representation of J' is 

J"{i) =^ Jni^^ mt jnt <! / (*(7), D) + h,{<^{-f))} U(d7) 



(4.13) 



7=(p,0)ex 



$(7)*(7)U(rf7) = £, / $(7)U(d7) = D 
7ex J7ex 



□ 



The proof of this is given in Section [S] The point of the second representation (|4.13p is that the innermost 
minimization problem (the one with $ and D fixed) involves linear constraints. Namely, that $ takes values 
in [0, 1] and that two integrals of $ take specific values. This will be useful in some of our numerical studies 
in the next section. 



5. Examples and Discussion 

In Subsection 15 . 1 1 we present some numerical examples. These examples showcase some of the possibilities 
and some of the implications of the dependence of the recovery distribution on the defaults. We conclude 
this section with Subsection 15.21 where we summarize our conclusions. 

5.1. Numerical Examples. Let's see what our calculations look like in some specific cases. To focus on 
the effects of recovery, let's assume a common (relatively high) probability of default of 8%; i.e., p^^" = 0.08 
for ah iV G N and n e {1, 2, . . . N}. 

In the first group of examples (Cases 1-4), we will consider four specific cases, one with fixed recovery rate, 
two homogeneous pools with stochastic recoveries and one heterogeneous pool with stochastic recovery. For 
comparison purposes, the law of large numbers average loss in the pool will be the same in all cases and equal 
to X ~ 0.064. We would like to remark the following. In a more realistic scenario one would expect to have 
pools with combinations of high-rated and low-rated counter-parties that have different default probabilities. 
Our formulation allows for such a scenario, but since we want to focus on the effects of recovery we assume a 
common probability of default in all of our examples. Moreover, recall that we have assumed as a reference 
point that all bonds have par value of $1. Wc will see that the tails (the large deviations principle rate 
functions) are significantly different. Although our theory has primarily focused on the rate function in the 
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large deviations principle for (Ljv)AfeN, the solution of the variational problem (|4.7p (or equivalently (|4.12p 
or (|4.13p ') gives useful information. 

In the second group of examples (Cases 5-6), we consider two heterogeneous pools, one with recovery 
whose distribution depends on the default rate (Case 5) and one with recovery whose distribution does not 
depend on the default rate (Case 6). As expected, the dependence on the default rate affects the tails. 

In Cases 2-5 we consider stochastic recovery rates. We want to consider the case that the recovery is in an 
appropriate sense negatively correlated with the defaults; i.e., that more defaults imply less recovery. In an 
economy that experiences recession, recovery rates tend to decrease just as defaults tend to increase. This 
property is also a documented empirical observation, e.g., see jSH09| . jABRSOS] and the references therein. 

In Case 1, let's assume that the recovery rate is fixed at 20%; i.e., g^'" = So,2 for all G N and 
71 e {1,2... A^}. The LDP here is essentially a straightforward application of Cramer's Theorem. Let's see 
how our general formulation of Theorem 14.61 covers this as a special case. In this case 7 = (0.08,^0.2) and 



U = Also we have 



and 

A^{e) ^ O.80iy{X) and A*(^) 



00 else 



if £ = 0.8i'(X) 
00 else 




Collecting things together, we have that (Z7v)ArGN and (ijAr)ArgN are governed, respectively, by the rate 
functions 

1 00 

We note that I[{£) is finite only if < ^ < 0.8. 

In Case 2 we consider a homogeneous pool with the recovery rate following a beta distribution. Essentially, 
this is a special case of Examples 13.11 and 14.41 For /3 > 0, define 

^^^{A)''^' 13 f {l~r)^-^dr, Ae^[0,l]. 

JreA 

this is the law of the beta distribution with parameters 1 and /?. As f3 increases, the amount of mass near 1 
decreases. We also have that 

/ rfip{dr) = — - 
Jre[o,i] t-hp 

for all /3 > (as /3 increases, the mean of fip decreases). This will allow a number of explicit formulae for 
the expected recovery (given the default rate). 

We here assume that the recovery rate has a beta distribution whose parameters depend linearly and 
monotonically on the empirical default rate. Namely, if the default rate is D, then the recoveries will all 
have common beta distribution with parameters 1 and 

/? = fMD) = 0.2 -0.1(£>- 0.08) ~ ^' 

note that / : [0, 1] ^ M+. Define 

g^siD, •) = Ai/.ff(D) for aU D G [0, 1] 

7aff = (0.08, gaff) 

Then U = (^-y^tf 

This choice of /aS results in a conditional expected recovery which is affine in D; i.e.. 



/ rg,ff{D,dr) =0.2~0.1{D-0.\ 

Jre[0,l] 



Observe that the expected recovery of a single name in the homogeneous pool is decreasing in D and always 
in [0, 1]. According to (|3.ip . the law of large numbers average loss is 

L = 0.08 X (1 -!- I = 0.08 X (1 - 0.2) = 0.064. 

V l + /aff(0.08); ^ ' 



The rate functions in this case are somewhat similar to those in Case 1. Again we have that 



oo . else 



For each /3 > 0, define 



Jr£[0,l] 



AUe) = sup {ei ^ Ap{e)} eeR. 

eeR 

Then \i v = j/(X)(5-y^j,. where v{X) > we have from Definition 14.31 

a:(£) = KX)A}^„.(,(,))(^) leR. 

Here {Ziq)^^^ and (LAr)jvgN are governed, respectively, by the rate functions. 

(5.1) ' ' \oo 

/^W = ^inf J;.o.osP) + M^^.,.) (A)} 



else 



In Case 3. we replace /aS of Case 2 with one that results in a conditional expected recovery which is 
quadratic in D] this allows us some insight into the effects of convexity in the conditional expected recovery. 
We set ^ 

^"^^^"^ ^ 0.2 - o.i(i:> - 0.08) - o.i(i:> - o.os)^ ~ ^' 

Again we have that /q : [0, 1] — > M+, and we set gq{D, •) = /iy^j-^)), 7q = (0.08, Pq), and have that U = S^^ 
Here we have that 

[ rQq{D, dr) = 0.2 - 0.1(1) - 0.08) - Q.l{D - 0.08)^ 

Observe that the expected recovery of a single name in the homogeneous pool is decreasing in D and always 
in [0,1]. 

We again get that L = 0.064. The corresponding rate function is /g. 

Case 4 involves the beta distribution again. Here, however, we now consider a heterogeneous pool of two 
types (Examples 13.21 and 14. 5p . We concentrate on the effect of the heterogeneity in the recovery distribution, 
so as in the previous cases all bonds will have default probability of 8%. For all D € [0, 1], every third bond 
will have recovery distribution governed by giaff and the remaining bonds will have recovery distribution 

average loss is: 



governed by g^. We thus have that U = + l^^^. It is easy to see that again the law of large numbers 



0.08 X i X I 1 5- -I + 0.08 X - X I 1 \ | = 0.064. 

3 V l + /aff(0.08)y' 3 V l + /q(0.08)/ 



For notational convenience and in order to illustrate the usage of Theorem 14.81 we use the alternative 
representation (|4.13p . If g{D,-) = for some / € C([0, 1]; M-|_), then for all D € [0,1], we have that 

Mg{9, D) ^ Af{D){S) for all 6* e R and Ig{£, D) ^ A*(^j(£) for all £ e [0, 1]. Since the support of U is exactly 
{7afi,7q}, we can consider $ and ^I* in B = -B(X; [0, 1]) of the form 

"5 = 4>AX{i^a} + <^SX{7q} and * = 1pAX{^.^ti] + '^BX{^^}- 
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Thus (|¥T^ becomes 

f 1 2 

4W = inf , inf inf <^ -D) + -^b/^JV-b, -D) 

0.4/3+20B/3=n 

1 2 

(5.2) 

f 1 - 2 - 

<t>Ai'A/3. + 2<t,ATpB/'i=e 

0.4/3+20B/3=n 

1 2 

In Figure we plot the rate functions /(, /j, and We use a Monte Carlo procedure to compute A 
and A*. As expected, all action functions are nonnegative and zero at the (common) law of large numbers 
average loss of Z = 0.064. We observe that < < I2 < I'l- In particular, the heterogeneous case, which 




Figure 1. Rate functions for fixed recovery, for the homogeneous cases and the heteroge- 
neous case. 

is a mixture of an affine conditional expected recovery and a quadratic conditional expected recovery, is in 
between the two homogeneous cases (cases 2 and 3). Of course, we should not be surprised that the rate 
function in Case 1 is larger than that in Cases 2 through 4, there are in general many more configurations 
which lead to a given overall loss rate. An interesting observation that seems to be suggested by the examples 
considered here is that more convexity results in smaller values for the rate function. 

The discussion from now on will be formal. However, we believe that it offers some useful insights in the 
effect of default rates on the average recovery of large pools. 

Another useful insight which we can numerically extract is the "preferred" way which losses stem from 
defaults versus recovery. For each £ G [0, 1], let D*(^) be the minimizeiQ in the expression (|5.ip for I2, or 
alternately the expression (|5.2p for I'^ . {£) can be interpreted as the "most likely " default rate in the pool 
given that Ln ^ £. We make the following assumption. 

Assumption 5.1. Assume that for each £ G [0, 1] the aforementioned minimizer D^,{£) exists and is unique. 

Numerically finding {£) in the examples considered here, verifies that Assumption 15.11 holds in these 
cases. For £ > L and 6 > 0, the Gibbs conditioning principle |DZ98[ Section 7.3] implies that we should have 



■^Casc 1 is of course degenerate this sense; for a given loss rate £, the default rate must be very close to £/0.8. 
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that 



(5.3) 



N 



lim Fn{\Dn - > S\Ln > ^} = 0. 



In other words, conditional on the pool suffering losses exceeding rate £, the default rate should converge 
to D^:{t}. Using this information, we can then say something about how the average recovery of the pool is 
related to the optimal default rate Z)h.(^). Motivated by p.ip we write that 

Average Loss=Default x (1-Recovery) ^ Recovery=l age 

Default 

to find an effective recovery rate in terms of the loss rate and the default rate. This formulation quantifies 
the fact that losses are due to both default and recovery. For atypically large losses in a large pool of credit 
assets, we should combine this with the Gibbs conditioning calculation of (|5.3p . Namely, let's define 

7^(^) = 1 - ^ 



[0, 1] [0, 1] be such that 



We can also formalize the dependence of recovery on default by letting TZ* 

(5.4) 7^* (£»*(£)) = n{i) for all £ e [0, 1]. 

We refer to 72,*(-) as the effective average recovery of the pool. 

Figure 5.2 is a plot of TZ* for the cases which we are studying. As it was expected, we observe that, for 
Cases 2, 3 and 4, the optimal average recovery of the pool and the optimal defaults are negatively correlated. 
Of course, this is consistent with the corresponding negative correlation of the individual recovery and 
default rate that is embedded in the choices of the individual recovery distribution. Notice however that TZ* 
is a "global" quantity in that it represents the effective average recovery in the pool. Another interesting 
observation is that the graph of the heterogeneous case which is a mixture of an affine and a quadratic 
conditional expected recovery is between the graph of the homogeneous cases 2 and 3 which treat the affine 
and quadratic case respectively. 




1 — Fixed Recovery 

2— Homogi 

3— Homogi 

4— Heterogeneous Affine-Ouadratic 



Figure 2. Effective average recovery of the pool versus optimal defaults. 



We conclude this subsection with a comparison of a heterogeneous portfolio whose distribution of recovery 
depends on the default rate with one that docs not. 

Case 5. Consider the setup of Case 4 with the only difference being the definitions of fa.ff{D) and fq{D). 
In particular, we consider 



fMD) 



1 



1, 



0.1 - 0.05(i:»- 0.08) 
•^''^^^ ~ 0.25 - 0.1(1? - 0.08) - 0.1(1? - 0.08)2 
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The corresponding rate function is Jg . 

Case 6. Again, consider the setup of Case 4 with the only difference being the definitions of fns{D) and 
fq(D). In particular, we consider 

/aff 

■'^ 0.25 

Notice that in this case the distribution of the recovery is independent of the default rate D. The corre- 
sponding rate function is /g. 

In both cases the law of large numbers average loss is the same as before. In Figure 5.3, we plot (a): on 
the left, the rate functions /g and /g and (b): on the right, TZ* for cases 5 and 6. 



Figure 3. Rate functions and TZ* when recovery distribution depends on defaults and when 
it docs not. 

As it is indicated by the left figure, one of the effects of the dependence of the distribution of the recovery 
on the defaults is to decrease the values of the rate function, i.e., Jg < /g. The effective average recovery, TZ* , 
is in both cases negatively correlated with the optimal default rate . Notice however that for large overall 
default rates the effective average recovery of the heterogeneous pool with dependence on defaults (Case 5) 
is less than the effective average recovery of the heterogeneous pool without dependence on defaults (Case 
6). This is consistent with our intuition, namely that dependence of the recoveries on defaults should affect 
the recovery of the pool and in particular that more defaults should decrease recovery. 

5.2. Conclusions and Discussion. In this subsection we summarize our findings and pose some question 
that would be interesting to study. 

• Assuming that the recoveries for the defaulted assets depend upon the number of defaults in a fairly 
general way we have characterized the typical (Lemma 13. 5p and atypical (Theorems 14.61 and 14. 8p 
behavior of the loss rates in the pool. This allows for consideration of the case when recovery rates 
arc affected by the number of defaults. We prove in Theorem 14.61 the large deviations principle for 
the joint family of random variables (Xat, vn) where Ln is the average loss and vn is the empirical 
measure on type-space determined by the names which default. Then, the large deviations principle 
for Ljv follows by the contraction principle. Furthermore, we derive in Theorem l4.8l various equivalent 
representations for the rate function which give some more insight to the favored way to rearrange 
recoveries and losses among the different types and also ease its numerical computation. 

• Moreover, as we have demonstrated in Subsection 15.11 the rate functions can be calculated numer- 
ically for both homogeneous and heterogeneous cases. The rate function determines the main term 
in the logarithmic asymptotics of the probability that the average loss in the pool is, let's say, bigger 
than a specific level. In particular, it determines the main term of the tail distribution of total losses 
on a portfolio consisting of many positions. See also |DDD041 IG02j for a related discussion. Also 
note that in the examples considered here, more convexity resulted in smaller rate functions. 
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• Our formulation allows to extract some information regarding the optimal default rate in the pool 
for a given level of average loss in the pool. Moreover, a useful insight into the fact that losses are 
due to both default and recovery is perhaps given by the effective average recovery as defined by 
relation ((O)) . 

• As it is indicated by the comparison of Cases 5 and 6, the effect of the dependence of the distribution 
of the recovery on the default rate is (a) : to reduce the rate function and (b) : to reduce the effective 
average recovery of the pool especially when the optimal default rate gets larger. 

Some interesting questions that naturally arise are below. These questions will be pursued elsewhere. 

• In this paper we deal with logarithmic asymptotics. It would be interesting to study the exact 
asymtpotics and characterize the prefactor that appears in the asymptotic expansion of the loss 
distribution as the pool gets larger. Similar questions have also been studied in jDDD041 ILMS09| 
under various assumptions. 

• A question that is in particular relevant for financial applications is to study measures such as 
Valuc-at-Risk (VaR) and expected shortfall (ES). VaR at level a G (0,1) is the a-quantile of the 
loss distribution and expected shortfall at level a is defined as ESa = E[L|L > VaRa{L)]. It would 
be interesting to study the asymptotic behavior of VaRaiLiq) and ESaiL^), characterize their 
limits as the pool gets larger and study the corresponding implications of the dependence of the 
recoveries on defaults. Under certain conditions similar questions have also been investigated in 
[G03llFM02llFM03] . 

6. Compactness of Level Sets 

The first part of the large deviations claim is that the level sets of / are compact. The proof follows along 
fairly standard lines. 

First, however, we need to topologize ./#i(X). This is done in the usual way. In particular, fix a point ★ 

that is not in X and define X+ =^ X U {*}. Give X+ the standard topology; open sets are those which are 
open subsets of X (with its original topology) or complements in X+ of closed subsets of X (again, in the 
original topology of X). Define a bijection t from ^i(X) to ,f3^(X"*') by setting 

{lv){A) v{A n X) + (1 - 6M) 

for all A S ^(X+ ) . The point * is introduced because is a subprobability measure. The topology of (X) 
is the puUback of the topology of ^(X+) and the metric on ^i(X) is that given by requiring t to be an 
isometry. 

Since X is Polish, so is X"'', and thus ,?3^(X+) is Polish, and thus ^i(X) is a Polish space. 
Proposition 6.1. For each s > 0, the set 

$(s) = {z e [0, 1] X .^i(X) : I{z) < s} 
is a compact subset of [0, 1] x ./#i(X). 

□ 

Proof. We first claim that $(s) is contained in a compact subset of [0,1] x .^i(X). Since [0,1] is already 

compact, it suffices to show that <I'a/(s) {i^ <= ^i(X) : H{h') < s} is a compact subset of ^i(X). If 
I' G ^m{s), then <C U and, since hp{x) = oo for a; > 1, we have that 

u{7eX:^(,)>l}^0, 

so for any A G ^(X), i'{A) < U(A). Since U itself is tight (it is a probability measure on a Polish space), 
<&m(s) is tight; for every e > 0, there is a CC X such that z^(X \ K^) < e for all i' e (f>j\/(s)- We claim 
that thus t(<f>A/(s)) is also tight. Indeed, fix e > 0. Letting io : X ^ X+ be the inclusion map, we have that 
io is continuous, and thus to(ife) is compact. Since singletons arc also compact, K* =^ Lo{Ki,) U {■*-} is a 
compact subset of X+. For every £ (f>jv/(s), ((i')(X+ \ A'*) = z/(X \ K^) < e, so indeed i($A/(s)) is tight. 
Thus $A/(s) CC ^(X+) and hence 

$m(s) C i-lt($A/(s)) CC .^i(X) 
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the last claim following since i is a homeomorphism. Gathering things together, we have that <J>(s) is indeed 
contained in a compact subset of [0, 1] x ^i(X). 

We now want to show that $(s) is closed, or equivalently, that ([0, 1] x ^i(X)) \ $(s) is open. Using 
Lemma HTTl wc have that 

([0,l]x^i(X))\$(s) 



<7ifEC(X) 



For each 9 e R and (p •= C'(X), then map (^, i') i-^ 91 + J^^y_4>{l)i^{d,j) — A^(6') is continuous, so we have 
written ([0, 1] x ./#i(X)) \ $(s) as a union of open sets. □ 

7. Large Deviations Lower Bound 

We next prove the large deviations lower bound. As with most large deviations lower bounds, the idea 
is to find a measure transformation under which the set of interest becomes "typical" . In this case, this 
measure transformation will come from a combination of Cramer's theorem and Sanov's theorem. 

First, we mention an auxiliary approximation result which will be useful in the proof. Its proof is in 
Section [TUl 

Lemma 7.1. Fix v S ^i(X) such that H(iy) < oo. Then there is a sequence (z^iv)neN such that 

lim = V and lim H{vn) ~ H{i') 

Af— >oo A— >oo 

[and thus vn ^ U for all € N) and such that 

di>N 



7 H> 



(7) and 7 = {p, g) ^ X(G,i){pWp (^^(7)^ 



are both well-defined and in C(X) for all A^ G N. 

We start with a simplified lower bound where the measure transformation in Cramer's theorem is fairly 
explicit. For each v € ^i(X), we make the usual definition |DZ98| [Appendix A] that 

dom a: =^ {£ e [0, 1] : < 00}; 

this will of course be an interval; ri dom A* will be the relative interior of dom A*. 

Proposition 7.2. Fix an open subset G of [0,1] x ./#i(X) and z = {£,1/) E G such that I{z) < 00 and 
£ (z ri domA*. Then 

(7.1) lim ^ln¥N{ZN eG}>-I{z). 

□ 

Proof. The proof will require a number of tools. For presentation purposes we split the proof in three steps. 
In Step 1, we prove several auxiliary results and are then used to define and analyze the measure change that 
is used to prove the initial lower bound. In Step 2, we prove that under the measure defined by (|7.7p in Step 
1, the recovery rates for the names that have defaulted are independent and that the default probabilities 
are independent as well. In Step 3, we put things together in order to prove the initial lower bound (|7.1|) . 
Step 1. Since £ G ri domA*, there is a S K such that 

(7.2) KiO)^£ and A:(^) = 0A;(0) - A,(0) 

(see |DZ98| [Appendix A]). Let's now fix a relaxation parameter j] > 0. Then there is an rji G (0,77) and an 
open neighborhood Oi of 1/ such that {£ — 7]i , £ + rji) x Oi C G. Using the first equality of (|7.2p . we have 
that {Al{9),i') = {£,!/) e {£- r]i,£ + 771) x d. Since the maps (77, £>) (A'i^iO) + ??, i>) and i> ^ Ao{9) are 
continuous, there is an 772 G (0, 1) and an open subset O2 of ^i(X) such that 

+V,i^)-fl& (0, 772), 7> e O2} C (£ - 771, £ + 771) X Oi 

^ ■ ' \A^{9) - A,{9)\ <77 fori>e02. 
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We next want to use Lemma 17.11 to choose a particularly nice element of O2 ■ Namely, Lemma 17.11 ensures 



that there is a z^* G O2 such that z^* <C U and such that both and 



(7.4) 



are in C(X) and such that \H{i'*) — H{i^)\ < -q. Let ©3 be an open subset of O2 which contains v* and such 
that 



7ex 



7ex 



< V 



for all !> e O3. 

We can now proceed with our measure change. For each N N, define 



N 



(7.5) ^ 9Lm (6) and ' ^ E A„0(7^^") " A,«,„ (^(7^-)) 

n— 1 

Standard calculations show that 



(7.6) 



E 



N 



cxp 



NA\ 



(N) 



= 1 and E 



N 



exp 



NA. 



{N) 



= 1. 



Define a new probability measure as 



(7.7) 



XACXp 



N ■ 



This will be the desired measure change. 
Define 



(7.8) 

On Sn, 



Sn = {\Ln - KJO)\ < 772, m e O3} . 
{Ln, vm) - {Kr. (^) + {Ln - {9)] , vn) e (£ - 7?i, ^ + 771) X Oi c G 

}] 



so in fact Sn C G. Thus 

Pjv{^w e G} > P^ 
Let's also assume that iV is large enough that 



XSn exp 



I 



7=(p,e)ex 



Ap(0(7))U(d7) 



Ap(0(7))Uw(d7) 

^7=(p,e)6X 

Thus by ([721) , (US]) , ([73]) and (H^ we have that 

4"^^ <e£ + \e\m - K{e) + ?; = + (|0| + 1) 

Next we prove that 



< V- 



(7.9) 



(7)) =^(7)0(7)-A,(0(7)) 



for U-almost all 7 = [p, g) e X. This follows from standard convex analysis and the form (|7.4p of (f> when 
p e (0, 1). Since H{v*) < 00, (|10.2p implies that, except on a U-negligible set. 



dv* 



(7) ) = fipip) = = p X - Ap(0) = ^(7)'/'(7) - Ap(0(7)) 



if 7 = (p, ^)) G X is such that p E {0, 1}. In other words, (|7.9p holds except on a U-negligible set. 
Therefore, 



(7.10) = / Hl)'^Nidj) 

J76X 



Ap(0(7))UAr(d7) 



7 = (p,e)6X 



< / 0(7)i^*(rf7)- 
J76X 



7=(p,e)ex 
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Ap((/)(7))U(d7)+277 



''7=(p,e)ex 



0(7)^(7) -Ap(</)(7))|u(d7) + 2?7 



A=(p,e)ex 



The last line of (|7.10p follows from 
Thus 

(7.11) Vm{Zn €G}> Vn{Sn) exp [~N {/(z) - {\0\ + 4) rj}] 

Step 2. Let's understand the law of {£n}i<n<N under PAr{-|^} defined by ()7.7p . In particular we prove 
that under the measure defined by ()7.7|) the recovery rates for the names that have defaulted are independent 
and that the default probabilities are independent as well. 

For any {^/'}i<n<Ar C M, we have that 



E 



exp 



N 
n=l 



N 

n 



A. 



re[0.1] 



exp 



Thus 
where 



exp 



AT "I "I JV [- 

^ ^ n ^ ^" / '^^P b/^M^ - r)] g7V,„(i'Jv(X), rfr) + 1 - A„ 

«=1 J J n=l[ -JreiO,!] 



for all iV G N and n £ {1,2... N}. In other words, the recovery rates for the names which have defaulted 
are independent with laws given by the gjv.„(i^Ar(X), •)'s. In particular, 

F \T 1^1 1 f . /,.,to,,(l-.)exp[^(l-r)]g^..„(.^(X),d.) 
[^-1^] -ivA.^" /„,,,,exp[^(l-.)]..,„K(X),d.) - 



Moreover, 
(7.12) Eat 



|iAr-A,„(0)nf^ 

1 " 



(1 - r)^eAr,„(t^Ar(X),dr) - / (1 - r)gAr_„(i/jv(X), dr) 



< 



1 

N' 



In a similar way, we next need to understand the statistics of the defaults under Pat. For {ijj}i<n<N C M, 

N UN 

= n {P""'" exp [V^i'n + 0(7"^'")] + 1 - P"'"} ■ 



E 



exp 



/^^{V^„A„ + 0(7^^")A„} 

ri=l 



Thus 

(7.13) 

where 



E 



•N 



exp 



TV 11^ 

n— 1 n— 1 



for all iV G N 
given by the p^'" 's 



and ?7 G {1,2... A'^}. In other words, under Pjv the defaults are independent with probabilities 
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Step 3. Let us go back to (|7.1ip . We want to show that hm y_ 
if hnijv-!.oo ^N{S'ff) = 0. To organize our thoughts, we write that 



N 



(Sn) > 0, which win in turn foUow 



(7.14) Fn{S;^)<PnWn^03} + En PN{\LN-KjO)\>f^2\&}x{.,eOs} 



V2 



E 



N 



Let us show that 
(7.15) 

Fix * e C(X). We write 
where 



1'(7)i/jv(d7) 



7ex 



*(7)i^*(d7) 



76X 



7ex 



'f(7)i^jv(d7) - / *(7)^*(d7) = £1 + £2 



76X 



1 ^ 

^=^l^{A„-p^-"}vI.(7^'.") 



ri=l 



n— 1 



v&(7)^*(d7) = / ^(7)vI/(^)U^.(d7)- / r_(^)vl,(^)U(d7) 



7ex 



7G 



X rfU 



7ex 



From (|7.13p we have that ]Ejv[^^i^] = 0; we also have by independence (Step 2) that 



E 



N 



\£i 



N\ 



< 



SUP7gX 1^(7)1' 

N 



The requirement that € C'(X) and Assumption 13.31 ensure that limjv^oof2 = 0. 

Combining things together, we have (|7.15p . Since ^1/ was an arbitrary element of C(X) and X is Polish, 
(|7.15p indeed implies (see |Str93| ) 



(7.16) 



N 



lim Fn{j^n ^ O3} = 0. 



Using (frT2|) and (frT6| we get by (fTTil) that hmw^oo Pjv(S'^) = 0. This and (fTTT]) give us the statement 
of the proposition. □ 

We can now prove the full lower bound 
Proposition 7.3. Let G be an open subset of [0, 1] x ^i(X). Then 

1 



lim — InPw {Zn G G} > - inf I{z) 

Ar-).oo ^ ^GG 



□ 



Proof. Fix z = {£,iy) e G. If I{z) < 00 and £ £ ri domA*, then we get (|7.ip from Proposition 17.21 If 
I{z) = 00, we of course again get (|7.ip . Finally, assume that £ £ domA* \ ri domA*. We use the fact 
that domA* C ri domA* and convexity oi £ 1-^ ^t{^)- Fix a relaxation parameter 77 > 0. Then there is an 
£' e ri domA* such that {£\u) e G and A*(^') < A*(£) + 77 (see [DZ98j [Appendix A]). Using Proposition 
17.21 we get that 



lim -IuPnIZn eG}>-I{£',,y)>-I{z)-ii 

Letting 77 \ 0, we again get (|7.ip . Letting z vary over G, we get the claim. 
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□ 



8. Large Deviations Upper Bound 



The heart of the upper bound is an exponential Chebychev inequahty. We will mimic, as much as possible, 
the proof of the upper bound of Cramer's theorem. The main result of this section is 

Proposition 8.1. Fix any closed subset F of [0, 1] x ^i(X). Then 



lim — IuPtv {Zn e F} < - inf Ilz). 

AT-J-oo TV z£F 



□ 



Not surprisingly, we will first prove the bound for F compact; we will then show enough exponential tightness 
to get to the full claim. 

Proposition 8.2. Fix any compact subset F of [0, 1] x .-#i(X). Then 



lim — InPjv {Zn e F} < - inf /(z). 



□ 



Proof. To begin, fix s < iniz^p Fix also a relaxation parameter 77 > 0. For each (0,0) G M x C(X), 

define the set 



= e [0,1] X .^i(x) : ee+ (^{^)u{d-i) - K{e) >s+ Xp{cp{^))u{d^) } 

(these open sets were used in the proof of Proposition 16 . 1 P . 

Fix now a z G F. By definition of / and Lemma HTTl we see that there is a {dz, fpz) G M x C(X) such that 
z e C(e,,0^). Since {£,1^) i-^ 0£ + J^^y_4'z{l)'^{dj) — A^{9) is continuous, there is an open neighborhood O* 
of z such that O* C Org and such that 



J + / Mi)Hdi) - MOz) >s+ Xp{Mi)Mdi) 



for all e O*. Thus 



zeF 



the compactness of F implies that we can extract a finite subset Z oi F such that 



z&Z 



and thus 



Fix now z & Z. We have that 



N 



{Zn e F} < ^ ^n{Zn e Ol}. 



z^Z 



N 



{Zn eOl}<PNi^0zLN + 



^z{l)vN{di) > s + K^iOz) + / Ap(</),(7))U(d7) 
7ex J7=(p,e)ex 



N 



exp [N {9zLn - {ez)}] exp 



N{ I Ul>N{di)~ / Ap(0.(7))U(d7) 
7ex Jj=(p^q)£x 



-Ns 

e exp 


H 











We have used here the fact that 



Xp{<PzmUN{dj) 



En [exp [N{ezLN - A,„(0,)}] \&] = 1 
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7=(p,e)ex 



Ap(0.(7))U(d7) 



(compare with (|7.6p ) and that 



E 



N 



exp 



7ex 



exp 



N / Ap(</>,(7))Ujv(d7) 



7 = (p,e)GX 



Letting — > oo, we get that 



lim ^lnPw{^jv e O*} < -s. 



This gives the claim. 

Let's next show that vj^ is in a compact set. 



□ 



Proposition 8.3 (Exponential Tightness). For each L > there is a compact subset ICl of ^i(X) such 
that 

Ih^^ ^IuFnIi^,, ^ ICl} < -L. 



□ 



Proof. First note that Assumption 13.31 implies that {UArjArgN is tight. Thus for each j e N, there is a 
compact subset Kj of X such that 



(L + j) 



n2 ■ 



We will define 



dcf 



ICl = <{ e ^i(X) : i^(X \ Kj) < ^— - for aU j eN 



1 

L + J 



Then ICl is compact, and we have that 
Wn ^ICl}<Y. 

Fn s vnO^ \ Kj) > — — 

oo 

< 5^ Pa, { iV(L + j^AT (X \ if, ) > iV(L + j) } 



< exp \-N{L + j)] Eat [exp [N{L + j)^mO^ \ K,)]] 



We now compute that 

Ea^ [exp [N{L + j)2z.A^(X \ A',)] ] = [] ^ 



JV 



TV 



exp 



N 



N{L + j)'-J2A^Xx\KAl''n 



exp 



N 



7V5]A,«,„((L + jrxx\K,(7^'")) 



< exp [iV(L + i)2UAr(X \ Kj)] < e 



AT 



We have used here the calculation that for 6* > 0, 

Xp{0)<ln{pe'> + {l-p)e'') = 
Combining things together, we get that 



JV 



j=i j=i j=i 



□ 



We can now get the full upper bound. 
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Proof of Proposition \8.1l Fix L > 0. Then 

Pn{Zn eF}< Pn{Zn eF,VNe JCl] + VNi^N ^ JCl}. 

We use Proposition l8.2l oii the first term together with the fact that [0, 1] x/C^ is compact. We use Proposition 
18.31 on the second term. Noting that (Lemma 1.2.15 in |DZ98| ) 

lim — In (Aaj + Sw) = maxj Um — InAiy, hm — InSiyl 

N^oc N N^oo N N^oo N 

we get 

li^ 4 InPwlZAr e i^, i^Tv e /ClI < - inf I{z) < - inf I{z). 

N^oo N z={l,v)<^F zeF 



□ 



9. Alternative Expression for the Rate Function 



In this section, we discuss the alternative expression for the rate function /' of Theorem 14.61 given by 
Theorem 14.81 In particular, this alternative representation shows that I'{t) has a natural interpretation 
as the favored way to rearrange recoveries and losses among the different types. In addition to providing 
intuitive insight, this alternative expression suggests numerical schemes for computing the rate function. We 
will rigorously verify that the alternative expression is correct, but will be heuristic in our discussion of the 
numerical schemes. 

We defer the proof of Theorem 14.81 to the end of this section and we first study the variational problem 
(|4.12p using a Lagrange multiplier approach. Even though an explicit expression is usually not available, 
one can use numerical optimization techniques to calculate the quantities involved. In order to do that, we 
firstly recall that we can rewrite J' of (|4.12p as a two-stage minimization problem, see expression (|4.13p . 

This naturally suggests an analysis via a Lagrangian. Define 



Ai, A2) - / {^{p)Is (vI/(7), D) + ftp($(7))} Wl) 



-AiM $(7)*(7)U(d7)-£!>-A2 
/76X 



/ cj>(7)U(d7)-i?| 



Let's assume that $* and ^* are the minimizers. Let's also assume that Ig{-,D) is differentiable for all 
7 ~ (Pi q) in the support of U. We should then have that for every r]i and 7^2 in ^ ~ B{X; [0, 1]), 

771(7) {/, (**(7), D) + ri;(<i>*(7)) - AiVl/*(7) - A2} U(d7) = 

7=(p,e)6X 

mhWh) (**(7),^) - Ai} u(d7) = 

7=(p,(?)ex 

Ignoring any complications which would arise on the set where $* = 0, we should then have that 

/, i^*h),D) + ^;($*(7)) = AiVl/*(7) + A2 
/;(vI'*(7),i?) = Ai 

for all 7 = {p, q) G X. This is a triangular system; the first equation depends on both Ai and A2, but the 
second depends only on Ai. Recalling now (j4.2[) and the structure of Legendre-Fenchel transforms, we should 
have that 

M;(Ai,I?) = vI;*(7) 
ft;(**(7)) = A2 + AiA/;(Ai,i?) - /,(**(7),Z?) = A2 + M,{\i,D) 
for all 7 = (p, q) € X. We can then invert this. This leads us to the following. Define 

dcf pe^2+A/,(Ai,_D) 

■^x,.D[Q) = Ki^^^D) Ai eM, C([0,1];^[0,1]) 
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where (Ai,A2) = {\i{£, D,{)), X2{i, D,(J)) is such that 

$Ai.A„Z5(7)U(d7) = ^ 



/76X 

We conclude this section with the rigorous proof of the alternate representation. 
Proof of Theorem \4.8\ First, we prove that J'{i) > !'{(-)■ Consider any $ and e S such that 

(9.1) / $(7)*(7)U(d7) = £. 

For any 6* G M, 

*(7)/, ( *(7), / *(7)U(d7) ) + hmi)) \ U(d7) 

7=(p,e)ex I V J7=(p,£,)GX / I 



$(7) sup J ^'^'(7) - A/e (e', [ 



$(7)U(d7) +^^($(7)) U(d7) 
7=(p,(?)ex [ e'eK \^ J^=(p,e)ex /J J 

> / 1 $(7) I 0f (7) - A/, ( / $(7)U(d7) ) I + ftp(*(7)) I U(d7) 

"'7=(p,e)ex I I V J7=(p,e)ex / I I 



$(7)A/, 0, / $(^)U(d7) ) + ;ip($(7)) \ U(d7). 
7=(p,e)ex L V J7ex 



Define v e .^i(X) as 



v{A) "^^^ / $(7)U(d7); A e ^(X) 
"'7eA 



then 



*(7), / '!>(7)U(d7) ) + M'^il)) } U(d7) 
7=(p,e)ex L V J7ex 



>9i- I M,ie, u{X))v{d^) + I h,[ ^(7) ) U(d7). 

'7=(p,(?)ex J7=(p,e)ex 



Varying 9, we get that 



$(7)/, *(7), / 'i>(7)U(d7) + /ip($(7)) U(d7) > + H(i^) > I\£) 

7=(p,e)ex 1^ \ J7=(p,e)ex / J 

and then varying $ and in i3 (such that (|9.ip holds), we get that J'{£) > I'{£). 
To show that !'{£) > J'{i), Gx u G ^i(X) such that < U. We want to show that 

(9.2) A:i£) + Hi,y)>j'i£). 

For all g € C([0, 1]; ^[0, 1]), define 

dof . 



a_{g, D) = inf{l -re [0, 1] : r G supp £>(£), •)} 
Dominated convergence implies that 



a+{Q, D) =' sup{l -re [0, !]:?'£ supp g(£), •)}. 



lim A;,(6I) =a„ =^ / a_(e, :/(X))i/(d7) 



7=(p,e)GX 



lim K'^{e) = a+ / «+(£», j/(X))i/(d7). 

From (j4.1ip and the monotonicity of moment generating functions, we can see that is nondccreasing; 
thus (a_,Q;+) e A'j,(M). This leads to three possible cases. 
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Case 1: Assume that i G (a_,a+), and let € R be such that A^(6'*) = £; i.e. 
(9.3) / M'^{9*,,y{X)Hd^)=i 



'7=(p,e)6X 

Then 



Alii) + H{y) = sup\9e- [ M,{9, v{X))v{d^) \ + f (^il)) U{dj) 

>9H- [ M,(9\ v{X))v(d^) + f hp (^h)) U(d7) 



dU 



7=(P:e)ex j'y=(p,g)eu 

/,(M;(r,KX)),KX))Kd7) + / f^(7) ) U(d7)- 



Define now $(7) =^ ^(7) and ^(7) M^(6'*, i^(X)). By (|T031) and (|4l0| respectively we have that 
$, * e S. Then dSS]) is exactly that / ^ $(7)*(7)U(d7) = £. Thus 



7=(p.e)GX 



$(7)/, (^«'(7),^^^$(7)U(d7)) +^^^($(7)) 



U(d7) > J'(^). 



This is exactly 

Case 2: Assume next that £ e [a+, 1]. For every g e C([0, 1]; i3^[0, 1]), define 

£1(9) '^^^ M,{9,i^{X)) - 9a+ig,iyiX)) = In f e-^("+(^■'^W)-(l-'^»^)(l^(X), dr). 
for all 9 eR; thus 

M,(^^,KX)) = 0a+(^,,l.(X))+£^(^^) c([o,i];^[o,i]) 



"'7=(p,e)ex 

for ah 6* e M. For ah e C([0,1]; ^[0,1]) and (1-r) e supp g(i^(X), •), the mapping 6* e-''("+(e,>'(x))-(i-r)) 
is decreasing and maps [0,oo) into (0, 1]. Monotone convergence implies that 

lim £U0) = lne{i^(X), 1 - a+{g, iy{X))}. 
0^00 

If ^ > a+, then we can use the fact that /^^(^ e)ex ^+(^)^('^'^) — ^ ^^-"^ all 6* > to see that 

A* (^) > \9{e-a+)- [ £l{9)v{d-i) \ > IW {9{e - a+)} = 00 > J'(£). 



If ^ = a+, then by the monotonicity of the f_^'s, 

(9.4) I,{a+{g,iy{X)),iy{X))=snp{-£l{9)}= Jh^ {-£l{9)}=ln 

eeWL e^co"- ^ \g{v{X),l - a+(g,iy(X))} 

for all g G C([0, 1]; ^[0, 1]), and 

Atil) = sup 1 - / f^(e)i'(d7) 1 - ^li^ I - / ft (e)i'(d7) I 

In f J 1 \ muV(^^)= / Ig{a^{g,v{X)))v{d^). 
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Defining $(7) =^ and ^'(7) = a+{Q, v(X)), we liave tliat 



<I>(7)vI/(7)U(d7) = / a+{Q. H^)Hdl) = e. 

7eX Jy = {p^g)l£X 



Collecting things together, we see that if ^ = a+, we again get 

Case 3: We finally assume that £ € [0,q;_]. The calculations are very similar to those of Case 2. For 
every g € C([0, 1]; ^[0, 1]), define 

Jr6[0,l] 

for aU 6* e M, so that 

M,{e, u{X)) = 9a_{g, i.(X)) + £^0)9 e C([0, 1]; ^[0, 1]) 

for all 6* e M. For aU g e C([0, 1]; ^[0, 1]) and (1 -r) € suppe(i^(X), •), the mapping 9 ^ e^«i-'')-"-(e''"W» 
is increasing and maps (—00,0] into (0, 1]. Monotone convergence implies that 

lim £^_{9) = \ng{i^{X),l-a-{g,iy{X))}. 
If ^ < Q;_, then we can use the fact that /^^(p e)ex'^-(^)^('^^) — ^ '^'-'^ all 6* < to see that 

A*(£) > li^ \e(e~a^)- [ £'LiO)u{dj)\ > li^ {6'(£-a_)} = 00 > J'(€). 

If £ = then by the monotonicity of the fi's, 



/,(a_(f,,i^(X)),i.(X)) = sup = lim {-Ei{e)] =\n(— ] 

foraU ge C([0, 1]; ^[0, 1]), and 



-a.{g,,y{X))} 



K{1) = sup - / em^idl) = lim - / Sm^id^) 

em I J7=(p,e)ex I 0^-00 1 J^^(p^g)gx 



1 

7=(p,e)ex'" V£'{i'(X), 1 - a^(g,v{X))} ) ' " J^^(p^g)(zx 



In ( „r...vN . . ■..v^^^ ) ^(^7) = / /,(a_(^., z.(X)))i.(d7). 



Defining $(7) =^ ^(7) and ^^(7) = «_ (g), zy(X)), we have that 



$(7)*(7)U(d7) = / a-{Q. H^)Hdl) = i 



again implying (|9.2 

Collecting things together, we have dHHl). We get ^^^^ by definng D = /^^^ * (7)^(^7)- Note that 



since $ and both take values in [0, 1], 

*(7)$(7)U(d7) < / $(7)U(d7). 
7ex J7ex 

This allows us to restrict the minimization in D to the interval [£, 1]. □ 
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10. Detailed structure of H 

In this section we prove Lemmas 17.11 and 14.71 The discussion in tliis section is somewhat technical. 
Recall the quantities defined in Definition 14.21 Note that 

(10.1) ho{x) ^ r " - " ° and hi{x) x ^ I 





Thus if H{v) < oo, we can restrict the region of integration to get that 



so in fact 



(10.2) 



7 = (p,e)GX 



U ^ 7 = b, f>) e X : p = and ^(7) ^0^=0 
U <i 7 = b, f?) e X : p = 1 and ^(7) ^l}=0. 



Fix ly G ./#i(X) such that H{i') < 00. The main technical challenges in both proofs is to stay away from 
the singularities in hp and h^. Note that 



and keeping (jlO.ip in mind, we thus need to be careful near p G {0, 1}, and for {x,p) G {0, 1} x (0, 1). 
To start, let's note some implications of the assumption that i?(i^) < 00. Clearly <C U. Secondly, 

(10.3) u(7GX: |;(7)>ll=0. 



dU 

Let's now do the following. Fix N ^ N. Define 

d^f J f (7) if P e 1 - ^] aud f (7) G (^, 1 - i) 



(10.4) C^(7) 



i ifpe [i,l-^] andf (7)<^ 

ifl^e [^,1-^] andf (7)>1 ' 



N 

Clearly < ^jv < 1, so we can define G ^i(X) as 

,yN{A) [ ^jv(7)U(d7)- A G ^(X) 

In light of (|10.3p and (|10.2p . limAr_>oo = ^ U-a.s., so it follows that limAr_^oo zy^r = i^- We next compute 
that 

hp (^(7)) if ^ < p < 1 - ^ and f (7) 
hp{j,) if ^<P<l-^andf (7)<^ 

hpil^i) if ^<P<l-^andf (7)>1-^ 
Using again (|10.3I) and (|10.2I) . we have that limjv-i-oo ^p(^Af (7)) = (lu^^)) ^'^^ U-almost-all 7 = (p, ^') G 
X. If p G [-^il";^]; then hp is increasing on [p,l] D [f ^"^'l] ^^^^ decreasing on [0,p] D [O' ■ 
Thus hp{^N{j)) < ^p(fu(7')) f^'" U-almost-all 7 ^ {p, g) G X. Dominated convergence thus implies that 
limjv->oo H{vn) = H{v). 

Proof of Lemma \7.1\ Fix N £ N; wc want to approximate by "nice" elements of C(X). Note that 

^^(7) = PX[0,1]\[N-\1-N-^{P) + X[N-\l-N-^]{p)^Nh)- 

Since U is regular (recall that X is Polish), we can approximate 7 = (p, g) 1— X[n-^.i-n-'^]{p)^n{j) by 
elements of C(X). From ([TUTi]) . we have that N^'^ < ^Nil) < 1 - if 7 = (p, £i) G X is such that 
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N ^ < p < 1 — N -'^j sowe can truncate these approximations at ^ and 1 — ^ without any loss. 
Namely, there is a sequence (C^'^)e>o hi C(X) such that 

^^^■^^ lim / \e''{l) ~ XlN-Ki~N-^]ip)^N{l)\ U(d7) = 0. 

For each e > 0, let ip, £ C([0, 1]; [0, 1]) be such that ip,{u) = 1 if u e [N-\ 1 - N-'^] and ip,{u) = if 
u e [0, 1] \ [N-'^ - e, 1 - N-^ + s]. For each e > 0, define 

for ah 7 = {p, g) E X. Then ^^.^ g (^(x) for all e > 0. We also have that 
/ |e''^(7)-^iv(7)|U(d7) 

"'7=(P,(?)6X 

< U{7 = e X :pe +£,1 +e] \ [N-\l~N-^]} 

+ I X[JV-i,i-w-i](p)|e'''(7)-Civ(7)|U(d7)- 
"'7=(p,e)ex 

Dominated convergence and (|10.5p then ensure that 



(10.6) hm / |r"(7)-ew(7)|U(d7)=0. 
Clearly < ^^'^ < 1, so we can define iyN,e € ^i(X) as 

Thanks to (|10.6p . we have that lim^^Q vn^^ = i^n. Note next that for 7 = (p, ^?) G X such that p € 

[0,l]\[iV-i-£,l-7V-i+e], 

f^pie'Hi)) - hpi^Nh)) = h,ip) - hp{p) = 0. 

If 7 = (p, g) € X is such that p € [N^^ - e, 1 - N^'^ + e], then 

(10.7) iV^l -£ < < 

so if £ < l/(2iV), 

|ftp(e''^(7)) - ftp(C^(7))| < ><k''^(7) -?iv(7)| 



where 

dcf 



H 



anp\\h!Jx)\ : < x < 1 - J— a.nd J— < p < 1 - J—\ . 



Thus if £ < 2^ , 



for all 7 = {p, g) e X. Thanks to (|10.6p . we thus have that limj_^o Hiy^,^) — H{v). 

We finally note that 7 ~ {p, g) i-> ^.^(^^'^(7)) is continuous on {7 = (p, gi) e X : p e (A^"-^ — £, 1 — 7V~-^+£)} 
f (|10.7p ensures that ^2 takes values in (0, 1) in this case). On {7 = (p, £>) G X : p e (0, 1) \ [N^^ — £, 1 — 
N-^ + £)}, we have that ?ip(^^'^(7)) = S^(p) = 0. This finishes the proof. □ 

Proof of Lemma \4.7\ Assume first that v is not absolutely continuous with respect to U. We will show that 
then the right-hand side of (|4.9p is infinite. Then there is an A e ^(X) such that i/{A) > and U(A) = 0. 
Since X is Polish, v is regular; i.e., 

iy{A) = sup {iy{F) : F C A, F closed} . 
Thus there is a closed subset of A such that z^(i^) > 0. Fix also now c > 0. For each n E N, define 

0ji (7) c exp [-n dist(7, F)] 7 e X 
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where dist(7,F) is the distance (in X) from x to F. Then c for all n E N, and 0„ \ cxf- Recall 

the definition of Xp{9) from (|4.8p . Since 6 i-^ Ap(6') is nondecreasing and continuous for each p G [0, 1], we 
also have that Xp{<pn{,j)) \ (0x^(7)) for all 7 = (j>, g) € X. Thus 

sup I / <t>{-i)u{dp)~ [ Ap(0(7))U(d7)l 



71— >00 



> \im { MlHdp)- / Ap(0„(7))U(d7) =ci.(F). 

'7ex •^7=(?',e)6X I 



Let c OD to see that the right-hand side of (|4.9p is infinite. 

Assume next that <C U. We use the fact that hp and Ap are convex duals of each other. For any 

^eC(X), 

HlMd-f) - / Ap(0(7))U(d7) 

To show the reverse inequality, let's write that 

H{u) = [ sup (0^(7) - \iO)} U(d7) - / lim F^(7)U(d7) 



where 



(10.8) ^^iv(7)= sup \ e—i^)~Xp{e) 



\e\<N 



for all e N and 7 = (p, gi) G X. We can explicitly solve this minimization problem; for N E N and 
7 = (p, e X, define 

(f (7)) if P e (0, 1), f (7) e (0, 1), and -N < h'^ (f (7)) < N 

(/)jv(7) = <(0 ifpe {0,1} and f5(7) =p 



dU > 

(where sgn is the standard signum function). Then 



^iVsgn(^(7)-p) else 



Fn{i) = 0^(7)^(7) - ^pi.4>N{.i)) 

for all 7 = (p, £>) G X and iV £ N. Clearly Fm and 0Ar are measurable, and (/)jv G B{X). From (|10.8p . we also 
see that Fn is nondecreasing in N. Thus by ()10.8p and monotone convergence 

H{iy)=lun f F^(7)U(d7) - lim / { Ml)^^) - \i<I^Nh))\ (^{dj) 



< sup M 0(7)^^(^7) - / Ap(0(7))U(d7) . 

ve-B(X) [J'y={p,e}ex J7=(p,e)ex J 

Since X is Polish, 1/ and U are regular; and thus we can approximate elements of B{X) by elements of C(X), 
completing the proof. □ 
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